Random Walks in Logarithmic and Power-Law Potentials, Nonuniversal Persistence, 
and Vortex Dynamics in the Two-Dimensional XY Model 
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The Langevin equation for a particle ('random walker') moving in d-dimensional space under an 
attractive central force and driven by a Gaussian white noise is considered for the case of a power-law 
force, F{r) ~ — r""^. The 'persistence probability', P(i{t), that the particle has not visited the origin 
up to time t, is calculated for a number of cases. For cr > 1, the force is asymptotically irrelevant 
(with respect to the noise), and the asymptotics of Poit) are those of a free random walker. For 
cr < 1, the noise is (dangerously) irrelevant and the asymptotics of Poit) can be extracted from a 
weak noise limit within a path-integral formalism employing the Onsager-Machlup functional. 

The case a = 1, corresponding to a logarithmic potential, is most interesting because the noise 
is exactly marginal. In this case, Po{t) decays as a power-law, Poit) ~ f~* with an exponent 6 
that depends continuously on the ratio of the strength of the potential to the strength of the noise. 
This case, with d = 2, is relevant to the annihilation dynamics of a vortex-antivortex pair in the 
two-dimensional XY model. Although the noise is multiplicative in the latter case, the relevant 
Langevin equation can be transformed to the standard form discussed in the first part of the paper. 
The mean annihilation time for a pair initially separated by r is given by t(r) ~ r^ ln(r/a) where a 
is a microscopic cut-off (the vortex core size). Implications for the nonequilibrium critical dynamics 
of the system are discussed and compared to numerical simulation results. 



I. INTRODUCTION 

This paper deals with two seemingly distinct topics: 
persistence and nonequilibrium critical dynamics. We 
shall show that they are related in the context of the 
nonequilibrium critical dynamics of the two-dimensional 
(2D) XY model with nonconserved order parameter. The 
link is established through the study of a d-dimensional 
random walker moving in a logarithmic potential. 

Persistence phenomena are related to first-passage 
problems for stochastic processes. Consider a stochas- 
tic process x{t). The 'persistence probability', Po(t), is 
the probability that x{t) has not crossed some given level, 
Xc (often taken to be zero), up to time t. The probability 
distribution, Pi{t), of the first-passage time (i.e. the first 
time t for which x{t) = Xc) is Pi = —dPo/dt. A familiar 
example is the ID random walk, with Langevin equa- 
tion dx/dt = ^(<), where ^(i) is a Gaussian white noise. 
For initial condition a:(0) = xq, the probability that the 
walker has not crossed the origin {x = 0) up to time t 
decays as Po(t) ^ x^lt^l"^ for i ^ Xq. The power enter- 
ing this algebraic decay has been termed the 'persistence 
exponent', Q, i.e. d = 1/2 for the ID random walk. 

Persistence phenomena have been widely studied in re- 
cent years p]-p^. Theoretical and computational stud- 
ies include spin systems in one |l[ and higher |2| dimen- 
sions, diffusion fields |^, fluctuating interfaces |4|, phase- 
ordering dynamics ||5[] , and reaction-diffusion systems [|6| . 
Experimental studies include the coarsening dynamics of 
breath figures M, soap froths ||], and twisted nematic 
liquid crystals PJ. Persistence in nonequilibrium critical 
phenomena has been studied in the context of the global 
order parameter, M{t), (e.g. the total magnetization of a 
ferromagnet), regarded as a stochastic process poi, but 



in the present work we will address a different, and more 
fundamental, aspect of nonequilibrium critical dynamics. 

If a system at its critical temperature evolves from a 
non-equilibrium initial state, critical correlations develop 
over a length scale, ^(i), which increases with time. Ac- 
cording to the standard theory [0, ^ ~ t^/^ for large 
t, where z is the critical exponent for equilibrium critical 
dynamics. This result, which is in accord with simple 
dynamical scaling, has been demonstrated within a field- 
theoretic framework flit]. This approach also shows that 
the result ^(i) ^ t^/^ is independent of the initial condi- 
tions. 

In the present work we challenge this picture for the 
specific case of the 2D XY model with nonconserved 
order parameter, and show that S,{t) ~ i^^^ if there 
are no free vortices present in the initial state, while 
^(t) r^ (t/lnt)^/^ if free vortices are present. Physi- 
cally, these two cases correspond to ordered initial states 
(e.g. the equilibrium state at T = 0), and disordered ini- 
tial states (e.g. the equilibrium state at T — oo) respec- 
tively. Furthermore, since the 2D XY model is described, 
through the Kosterlitz-Thouless (KT) theory ||l||, by a 
critical line, T < Tkt, rather than a single critical point, 
the above dependence on initial conditions will persist 
throughout the KT phase. 

The connection with persistence arises as follows. For 
an initial condition containing free vortices and antivor- 
tices, vortex-antivortex annihilation is the slowest relax- 
ation process. Much can be learned by studying the an- 
nihilation of a single, initially widely separated, vortex- 
antivortex pair. The probability that they have not an- 
nihilated up to time t defines a persistence problem. By 
a series of transformations, this can be mapped onto 
a random walk in a logarithmic potential. Analysis of 



this problem shows that the persistence exponent 6* is a 
continuous function of the ratio of the strength of the 
potential to the strength of the thermal noise, i.e. 9 is 
nonuniversal. In the context of the vortex problem, we 
find 9 = 7rps{T)/kBT, where Ps{T) is the large-scale 
spin-wave stiffness at temperature T. Since the ratio 
Ps/kBT takes the universal value 2/7r at Tkt, we obtain 
9{Tkt) = 2, while 9{T) ^ oo as T ^ 0. 

Although the vortex problem in two-dimensions was 
the initial motivation for the present study, for pedagog- 
ical reasons we will begin (section II) by discussing the 
d-dimcnsional random walk in a logarithmic potential, 
and deriving the persistence properties. We show that 
the exponent 9 depends continuously on the ratio of the 
strength of the potential to the strength of the noise [Q , 
i.e. the logarithmic potential is marginalin the Renormal- 
ization Group sense. In section III we consider the case of 
a general attractive power law force, F(r) = {—A/r°')r, 
where the logarithmic potential corresponds to ct = 1. 
We show that for ct > 1 the force is irrelevant (relative 
to the noise) as far as asymptotic persistence properties 
are concerned, so the results of the free random walk still 
hold. For cr < 1, by contrast, the force is relevant and the 
noise becomes irrelevant. For zero noise, however, the dy- 
namics are deterministic, so the persistence asymptotics 
are governed by rare fluctuations of the noise. In this 
sense, the noise is dangerously irrelevant. The asymptotic 
persistence follows from an optimal path (or steepest de- 
scent) approach formally valid in the weak noise limit 
Q,|5|. The result is the 'stretched exponential' decay, 
Po{t) ~ exp[-const.i(i-'^)/(i+'^)]. 

Section IV deals with the application of the previous 
results to the problem of vortex-antivortex annihilation 
in the KT phase of the 2D XY model. The potential en- 
ergy of a free vortex-antivortex pair is logarithmic in their 
separation, so this is the marginal case with a nonuniver- 
sal 9. An additional complication is that the effective 
vortex mobility is scale (i.e. separation) dependent, im- 
plying, via the Einstein relation between mobility and 
diffusion constant, that the noise strength is also scale de- 
pendent or, equivalently, that the noise is multiplicative. 
However, this system can be transformed to an additive 
noise problem identical to that discussed in section II. 
Although the persistence exponent is non-universal, the 
dependence of the typical annihilation time on the ini- 
tial separation r has the universal form t ^ r^ln(r/a), 
where a is a (non-universal) short-distance cut-off, for all 
T < Tkt- Standard scaling arguments then suggest that, 
for an initial condition containing many free vortices and 
antivorticcs, the characteristic length scale describing the 
approach to equilibrium will grow as ^(i) ~ (i/lnt)^" 
throughout the KT phase. 



II. RANDOM WALKER IN A LOGARITHMIC 
POTENTIAL 

The Langevin equation for an isotropic, continuous- 
time random walker moving in d-dimensions in the cen- 
tral potential V{r) = A\nr is 



dr A ^ ^, , 

at r 



(1) 



where the Gaussian white noise has the correlation func- 
tion (i,j = l,...,d) 

{^^mJ{t'))=2D^,,S{t-t'). (2) 



A. Reduction to a One-Dimensional Problem 

The analysis starts from the equivalent Fokker-Planck 
equation for the probability distribution, P(r, i), for the 
position of the particle at time t: 



dt 



DVP] 



This equation can be reduced to an effectively one- 
dimensional equation in terms of the radial coordinate 
r by integrating over the angle variables. Writing r = 
(r, i7), and defining the radial probability distribution 



Q{r,t) 



-^ J dnP{r,n,t), gives 



dQ 
dt 



Ddr -Q + drQ 



where dr = d / dr and 



h = 



A 
D 



1-d. 



This is equivalent to the radial Langevin equation 
dr _ A {d-l)D 
dt r r 



-m, 



(4) 



(5) 



(6) 



where {^{t)^{,t')) = 2D5{t-t'), i.e. projecting (|) onto 
the radial direction leads to an additional repulsive force 
{d—l)D/r, proportional to the noise strength. This mean 
[see (o)] that a particle moving in a logarithmic poten- 
tial in space dimension d has the same radial distribution 
function as a free particle (random walker) moving in di- 
mension d' = d — A/D (which need not be an integer, or 
even positive). 



B. The Persistence Probability 

In this subsection we derive an exact expression for the 
persistence probability Po(^)- First we make the change 
of variable Q = r^-^~'''^^^R in (Q). Then the equation gov- 
erning the relaxation modes R{r, t) — Rk{r) exp(— Dfc^t) 
becomes 



d'^Rk , IdRk 



dr"^ 



dr 



Rk = 0, 



(7) 



where 



v = {l + b)/2 . 



(8) 



The solutions are the Bessel functions Ju{kr) and 
J-^{kr), so the general solution is 

/>OC 

Q(r, t) = r(i-'')/2 / dk [a{k)J^^{kr) 
Jo 

+ /3{k)J^{kr)] exp{-Dk^t). (9) 

This desired solution has to satisfy the initial condi- 
tion (5(r, 0) = d{r — To). To determine Po(t) we impose 
an absorbing boundary at r = 0, such that the particle 
is removed if it reaches the origin. To determine the ap- 
propriate solution we note from (Q) that the probability 
current is 



j{r)^-D{drQ + bQ/r). 



(10) 



Consider separately the two terms in the general solution 
(y). The term involving J_j^ behaves as r for r — > 0, while 
the term in J^, behaves as r~^. Both terms have correc- 
tions which multiply the leading terms by power series 
in r^. Inserting both forms into ( pl]| ) shows that the first 
term gives a finite (and negative) current at r = 0, so 
that for this solution the origin is an absorbing point, or 
current sink. The second term, on the other hand, gives 
zero current and the origin is a not a special point: This 
solution therefore gives the spherically averaged Greens 
function for diffusion in d' = 1 — & dimensions. 

Using the orthogonality properties of Bessel functions 
to determine the functions a{k) and (3{k) in each case, 
we find that the Greens function is given by 



QG{r,t) =ro 



ro 



(l-b)/2 



dkkJ-i, (fcro) J-,y (fcr) 
X exp(-L»fc2t) 



2Dt I ro 



(l-b)/2 



exp 



ADt 



--"(S)' <"' 



where I,^{z) is the modified Bessel function, while the 
relevant solution for an absorbing point at the origin is 



]absir,t) = ro — 



(l-b)/2 »oo 



dk kJi, {krQ)J^ [kr) 
X cyi^{-Dk'^t) 



2Dt \ro 



(l-b)/2 



exp 



„2 I ^2 



ADt 
frrp \ 



At this point a comment on the possible values of b is 
needed. The leading term in Qg for r ^ is 0{r~''), 
and this gives rise to zero current at r = as discussed 



above. The next-to-leading term, however, is 0{r'^~^) 
leading to a current of order r^~''. This vanishes as 
r — > only for 5 < 1. Thus the Greens function is ill 
defined for 6 > 1: The non-zero current at r = in this 
regime means the radial probability distribution Qg (r, t) 
collapses onto the origin. In terms of the dimension d' of 
a free random walk we have b = 1 — d', so the require- 
ment 6 < 1 means d' > 0, which makes physical sense. 
A physically reasonable Greens function can be restored 
for the case d > 0, A > dD, where 6 > 1, by regulating 
the r — singularity of the force, F{r) = —A/r. For ex- 
ample, if _F = —A/(r + e) we expect the Greens function 
to have a width which vanishes with e. 

We turn now to the more interesting case, for present 
purposes, of the distribution Qabsir^t), appropriate to 
an absorbing point at r = 0. Equation (n3) gives 
Qabs ^ r ioY r ^ Q, so the current (|l(]| ) at the origin 
is :;'(0) = —D{1 + b){Qabs/r)r=o- This requires b > —1, 
since j(0) < for an absorbing boundary, and Q is nec- 
essarily non-negative. In the terms of the equivalent free 
random walk with dimension d', the condition fe > — 1 
requires d' < 2, i.e. the probability to reach r = for a 
random walk in d' > 2 dimensions is zero. 

For 6 > — 1, the persistence probability, Po{t), can be 
readily calculated from (|l^) by first computing the cur- 
rent (^0|) at r = 0. This current gives the rate of change 
of the persistence. 



^■<°)^^ 



-Flit) 



(13) 



where Pi (t) is the probability distribution for the time at 
which the particle first reaches the origin ('first passage 
time'). Using (O) and (p^ gives the final result: 



Pi{t) 



AD 



r[(l + fe)/2] r2 [Am J 



2 \ (3+b)/2 






exp 



ADt 



(14) 



The persistence probability is Po (t) = Jf ds Pi (s) . Us- 
ing the large-t behavior of ([l4| ) gives 



Po{t) 



1 



.^yi.b)/. 



r[(3 + 6)/2] \ADtJ 



(15) 



for t — > CX3. Thus the 'persistence exponent', 9, is given 

by 



(l + 6)/2. 



(16) 



and is non-universal. For the free random walk in dimen- 
sion d', this translates to = (2 — d')/2. 

The nonuniversality of 9 with respect to the strength, 
A, of the potential is special to the case of a logarithmic 
potential, for which the Langevin equation (m) is invari- 
ant under the rescalings r -^ ar, t — > a^t of space and 
time. This means that the potential is a marginal per- 
turbation with respect to the equation with yl = 0, and 



it is this marginality which is responsible for the contin- 
uous variation of with A (actually, with the ratio A/D) 
through it's dependence on 6, which we recall is defined 
by (|^). The condition b > —1 for the particle to visit the 
origin with probability one is equivalent to ^ > (d—2)D. 
Note that for d = \ this even allows a (sufficiently weak) 
repulsive potential, whereas for d > 2 a guaranteed visit 
to the origin requires a sufficiently strong attractive po- 
tential. 

In the following section we discuss the case of a gen- 
eral power law potential, corresponding to a force F{r) = 
—Ajr". We show that the force is irrelevant for cr > 1, 
and the asymptotic persistence probability is that of a 
free random walker. For cr < 1 the force is a relevant 
perturbation to the free random walker. In this case the 
noise term is irrelevant, but 'dangerously irrelevant' as 
far as the calculation of Paif) is concerned. We show 
that in this case -Po(i) decays as a stretched exponential. 



t = ro+'^/[(l + a)A\. The limiting behavior of Po(i) at 
large t is therefore dominated by rare events: the noise is 
a dangerously irrelevant variable in this context, and we 
cannot simply set it to zero. Instead, we have to examine 
the limit of small but non-zero D. 

We can argue for the asymptotic form of Pq (t) as fol- 
lows. A long survival time of the particle is a rare event, 
dominated by an activated process where the particle ini- 
tially moves away from the origin. The potential corre- 
sponding to the force —A/r'^ is V{r) = Ar^~'^ /{I — a). 
Suppose the particle is driven (by the noise) to a point 
ri. The time for a subsequent deterministic descent to 



the origin is t 



K.1+0- 



/[A(l -(- cr)]. The activation barrier 



for the 'uphill' process is AF = [Al{l — a)]{7 
so the probability to reach ri before the origin is of order 
exp(— [A/(l — a)]r\~'' /D)^ where the term in ro in AT^ 
has been taken out and absorbed into a pre-exponential 
factor. Using the time for the subsequent free descent to 
estimate ri, i.e. n ~ [A{1 + cr)t]^/(^+'^\ gives 



III. RANDOM WALKER IN A POWER-LAW 
POTENTIAL 

After reducing the problem, as in the previous sec- 
tion, to an effectively one-dimensional problem for the 
distance, r, of the particle from the origin, the radial 
Langevin equation reads 



dr 



A 



{d-l)D 



-m, 



(17) 



A. Scaling Analysis 

Under the scale transformations r —f ar, t — > a^t, the 
Langevin equation (O) retains the same form, but with 
potential strength A and noise strength D rescaled to 



D' = a'-^D. 



'A, 



(18) 
(19) 



For a =/= 1 these equations have two non-trivial fixed 

points: 

(i) A > 0, D = 0, with z = 1 + cr, which is stable {D 

scales to zero at large time) for cr < 1, and 

(ii) A — 0, D > 0, with z — 2, which is stable for cr > 1. 

For the special case cr — 1, there is a line of fixed points 

with z — 2 and A/D arbitrary. In this last case, as we 

have seen, the exponent 9 depends continuously on A/D. 

At the second fixed point, with cr > 1, the force is ir- 
relevant: it falls off too rapidly with distance to affect 
the asymptotic large-time behavior. In particular, the 
exponent 9 is given by the zero-force value, 9 = {2~ d)/2 
(0 < d < 2). 

At the first fixed point, with ct < 1, the noise is irrele- 
vant. If we set D to zero, the process (O) is determinis- 
tic: A particle staring at ro reaches the origin in a time 



Po{t) --exp 



A 



{l-a)D 



[A(l + a)i](i^'^)/(i+'^M , (20) 



a 'stretched exponential' form, Po(0 ~ exp(— const. t'^), 
with (3 — (1 — cr)/(l-fcr). Note that the coefficient of 
i'^ does not depend on the initial displacement rg. The 
reason is that the rare trajectories for which the particle 
survives a long time take the particle far from its original 
position. 

Despite its crudeness, this argument gives the correct 
result up to a constant of order unity in the exponent. 
This is shown in the next subsection using a path-integral 
formalism, augmented by a steepest descent calculation 
valid for Z) ^ 0. Finally we note that the crucial fac- 
tor (Ai)(i-'')/(i+'^V^ in tlie exponent of (|2^), and in 
particular the value of the exponent /3, can be deduced 
immediately from dimensional analysis once one recog- 
nizes that (i) the result must be independent of tq, and 
(ii) the factor \/D is a necessary consequence of activated 
dynamics. 



B. Path Integral Formulation for a < 1 

We begin from the probability distribution functional 
for the noise history, ^(i). Since the noise is Gaussian 
and white, this functional is 



P[^(i)]=^exp(^-^y'diC'w), 



(21) 



where TV is a normalization constant. This can be trans- 
formed to a probability distribution functional for r{t) 
using the Langevin equation (O): 



P[r[t)] =NJ[r{t)\ exp 



D 



S[r{t)] 



(22) 



where J[r] is the Jacobian of the transformation from 
^(t) to r{t), whose precise form will not concern us, and 
S[r] is the Onsager-Machlup functional (or 'action'), 



SM^/ 






(23) 



It is convenient to compute Pi (i) , the probability den- 
sity for the first visit to the origin (recall that Pi = 
—dPo/dt), given that the particle starts from tq at t = 0. 
This is given by the path integral 



Piit) 



dr{t) J[r{t)] e^p [ ~-S[rit)] 



(24) 



where the time integral in ( p3| ) now runs from to t 
(and we introduce a dummy time integration variable s), 
and the path integral is over all paths r(s) which satisfy 
r(0) = ro and r(i) = 0. The use of ~ in ( |24| ) means we 
are concerned only with the leading exponential terms 
and not with the prefactors. 

In the limit Z? — > 0, the path integral can be evalu- 
ated by steepest descents. To leading order the Jaco- 
bian, J[r{t)], can be replaced by J[rc{t)], where rc{t) is 
the 'classical' path which minimizes the action S[r]. The 
Jacobian therefore contributes to the prefactor, and we 
will not consider it further. In a similar way, the term 
{d — l)D/r in S is subdominant for D ^ and can be 
dropped to leading order. 

As noted above, the dominant path (path of least ac- 
tion) is the one which keeps ^{t) as small as possible for 
as long as possible. This is achieved by having the parti- 
cle move initially away from the origin ('uphill'), where a 
smaller noise force is needed to overcome the determinis- 
tic force driving the particle towards the origin. The path 
rc{t) therefore consists of two parts: an 'uphill' part to 
a maximum displacement ri, followed by a determinis- 
tic (^ = 0) downhill part. Only the uphill path has a 
nonzero action associated with it. 

The variational problem for the uphill path is sim- 
plified [ |lq | by introducing the velocity v = dr/ds and 
parametrizing the path by w(r) instead of r{s). The ac- 
tion then becomes 



^. . 1 n dr f A 
Sv]^- / — [v+ — 



4./ro « 



(25) 



The variational equation SS/Sv — becomes 1 — 
A'^/ivr")'^ = 0, with solutions v = ±A/r''. The mi- 
nus sign corresponds to the deterministic downhill path, 
with zero action, the plus sign to the uphill path. The 
action for the latter is 



S = AI ^ = 



A 



1 



(26) 



The times tu-, td associated with the uphill and downhill 
paths respectively, are 



r^ dr 

tu= — = 



1 



V A{1 + a) ' 



nl+cr 



td 



A{l + ay 



(27) 
(28) 



The final step is to eliminate ri in (|2^) in favor of the 
total time t = t^ + td- For t —^ cxd, ri ^ oo and tq <^ ri 
in (p7|). Dropping the terms in rp in (GO) and ( p7| ) gives 
n ~ [A(l+cr)t/2]i/(i+'") and 



Pi{t)^cxp{-S/D) 



exp 



A 



{l-a)D 



A{1 



(1-<t)/(1+<t)^ 



(29) 



The persistence probability, Po{t) = j^ dtPi{t), clearly 
has the same asymptotic form. It differs from (Efl) only 
by factors of order unity, as promised. It is easy to show 
that the corrections due to keeping the rg term in ([27) 
vanish for f — > oo, while the correction associated with 
the ro term in (Eq) represents a time-independent pref- 
actor. This prefactor, exp[Arl~'^ / {1 — a)D] can of course 
be very large (and very sensitive to rg) for small D. 

It is important to note that while the asymptotic form 
( p9| ) was derived in the limit 13 ^ 0, it actually holds as 
an asymptotic result for all D, since D is an irrelevant 
variable and scales to zero as i — > oo. 



IV. VORTEX- ANTIVORTEX ANNIHILATION IN 
THE 2D XY MODEL 

A. Nonequilibrium Critical Dynamics 

As a final application of these methods, we consider 
vortex dynamics in the two-dimensional (2D) XY model, 
which was in fact the motivation for the present study. At 
all temperatures T < Txti where Txt is the Kosterlitz- 
Thouless (KT) transition temperature, a system pre- 
pared in a nonequilibrium initial state will approach 
the equilibrium state through a coarsening mechanism 
in which local equilibrium is established over a length 
scale S,{t) which grows with time. For example, the spin- 
spin correlation function has, according to the KT the- 
ory iQ, the equilibrium form C{r) ~ r~'^^'^\ for all 
T < Tkt- Consider now a system prepared in a random 
initial state, with only short-ranged spatial correlations, 
and allowed to evolve in contact with a heat bath at tem- 
perature T < TxT- According to the conventional the- 
ory of nonequilibrium critical phenomena pT[ | , the system 
will approach equilibrium via a dynamical scaling state, 
characterized by a growing length scale ^(i). For exam- 
ple, the scaling form for the spin-spin correlation function 
reads 



C(r,i) 



J, J) 



/ 



m) 



(30) 



The scaling function f{x) has the Hmiting behavior 
/(O) = const, (so that equihbrium is recovered for £^{t) -^ 
oo), while f{x) falls off rapidly for x S> 1, representing 
the fact that the spins are uncorrelated on length scales 
large compared to ^(i). 

The standard theory of nonequilibriuni critical dynam- 
ics predicts that the length scale £,{t) should grow as 
^(t) ~ t^'^, where z is the equilibrium dynamical expo- 
nent. This should hold independent of the initial condi- 
tions, though the scaling function f{x) in ( pO| ) can depend 
on initial conditions. A commonly considered case is uni- 
form initial conditions (all spins parallel). For this case 
one requires that f{x) ~ x'' for a; — > oo, since the long- 
range order present in the initial condition will persist at 
any finite time. 

In a recent paper, however, Bray et al. Pq | have ar- 
gued, on the basis of numerical simulations and physical 
arguments, that this picture breaks down for the 2D XY 
model with nonconserved dynamics ('Model A' of the 
Hohenberg-Halperin m% classification). Specifically, the 
growth of ^(i) will depend on whether the initial state 
contains free vortices and antivortices. In particular, they 
argue that for a uniform initial condition, for which there 
are no free vortices, ^(i) is determined by spin wave the- 
ory to be ^{t) ~ i^/^ [ESJ. On the other hand, for a ran- 
dom initial condition there are many free vortices and 
antivortices present. The dominant coarsening mecha- 
nism in this case is vortex-antivortex annihilation, and 
this leads to £,{t) - {t/\nty/^. 



B. Vortex- Antivortex Annihilation 

Physical arguments for ^(t) ~ (t/lni)^/^ have been 
given previously for the coarsening dynamics at T = 
from a random initial condition [ p9pc| ]. The basic idea 
p9[ is to consider a single vortex-antivortex pair, and to 
derive expressions for the energy, Ep{r) of the pair, and 
hence the force, F{r) = —dEp/dr, between them, as a 
function of their separation r. The result is Ep '^ ln(r/a) 
(where a, the vortex core scale, is a microscopic length), 
and hence F ~ — 1/r. Crucially it is found |19] that the 
vortex mobility fi, which relates the force, F, to the ve- 
locity V, via V — fiF, also depends logarithmically on the 
pair separation: fj, ~ l/ln(r/a). At T = 0, therefore, the 
variable r obeys the deterministic equation 



In 



dr 
'di 



(31) 



up to an overall constant which can be absorbed into 
the timescale. Integrating this equation gives the anni- 
hilation time for a pair initially separated by a distance 
ro ^ a, namely i ~ Tq ln(ro/a), which one can invert to 
obtain tq ~ [i/ln(t/a^)]^/^. This already suggests that, 
in a many vortex situation, lengths and time are related 
by^(i)~(i/lni)'/'- 



This result can be motivated in another way using scal- 
ing arguments. If the characteristic scale in a many- 
vortex system is ^, the typical force on a vortex scales 
as 1/f, the typical mobility as l/ln(^/a) and the typical 
velocity as d^/dt. This gives ln(^/a) d^/dt ^ 1/^, i.e. 
^- (i/lni)^/^ as before. 

We now discuss the influence of thermal fluctuations 
on the annihilation of a single vortex-antivortex pair at 
non-zero temperatures. 



C. Vortex- Antivortex Annihilation at T > 

As a first step, we present a detailed and quantitative 
treatment of the T = arguments employed in the pre- 
vious subsection. It is convenient to adopt a continuum 
approach based on the non-linear sigma model Hamilto- 
nian 



H 



d^r{Vcf)f 



(32) 



where is the two-component order parameter field, sub- 
ject to a local constraint (fp' — 1^ and ps is the spin- wave 
stiffness. For a field configuration describing a single free 
vortex, (f) — r/H, one has (V(/))^ — l/r^, leading to an 
energy E^ = {ps/2) J{d^r/r^) = Trps\n{L/a), where L 
and a are the system size and microscopic cut-off respec- 
tively. A vortex-antivortex pair, separated by distance 
r, screen each other's far fields at scales larger than r, 
leading to a pair energy Ep(r) ~ 27rps ln(r/a), and an 
attractive force F = ~ dEp/dr — — 27rps/r between the 
vortex and the antivortex. 

The corresponding continuum description of the non- 
conserved ('model A') dynamics is given, at T = 0, by 
the Langevin equation [|7| 



JH 



dt 



(33) 



This equation can be used [|19| to compute the effective 
friction constant 7(r) = l//i(r), where p, is the mobil- 
ity, associated with the motion of the vortex and an- 
tivortex under the force F . An isolated vortex mov- 
ing at speed v in the x-direction has field configuration 
(f){x, y, t) — (j)fj{x — vt, y). Energy is dissipated at a rate 




-IvV 



(34) 



where the notation 7^ indicates the residual depen- 
dence of 7 on V at non-infinitesimal velocities. Inserting 



the equilibrium vortex configuration, which is isotropic, 
gives the hmiting zero- velocity friction constant: 70 = 
Ey/psT ~ (tt/F) ln(i/a), i.e. 70, like the vortex energy 
Ey, diverges logarithmically with the system size, L. For 
a vortex-antivortex pair, this translates into a logarith- 
mic dependence on the separation |19[| , 



7(r) 



f)'< 



(35) 



The effect of thermal fluctuations, neglected up to now, 
is twofold. Firstly, as in the equilibrium theory, thermally 
activated bound vortex-antivortex pairs lead to a renor- 
malization of the spin- wave stiffness, ps, and kinetic co- 
efficient, F, to temperature-dependent functions, Ps{T) 
and r(T), that describe the large length-scale proper- 
ties of the system. In equilibrium, however, there are no 
free vortices at any temperature below the KT transi- 
tion temperature, Tkt- This means that the large-scale 
properties are described by the spin-wave theory. In this 
theory all vortices are neglected and the angle represen- 
tation, (/)(r) = [cos0(r),sin0(r)] is employed, with the 
angles 9{y) defined on the interval (— cx),cx)). The ef- 
fective Hamiltonian for the long-wavelength degrees of 
freedom is H = [p,{T)/2] J (fr {V9)^, and the equa- 
tion of motion is 80 /dt = -r{T){6H/6e) + ^(r,t) = 
r(T)ps{T) V^6'-f ^, where ^(r, t) is a Gaussian white noise 
with correlator given by the fiuctuation-dissipation the- 
orem, (e(r, i)^(r', t')) = 2F(T) keTSir - r')S{t - t'). 

Since the equation of motion is linear, it can be solved 
exactly. The dynamic exponent is z = 2 (for all T < 
Tkt)- Nonequilibrium properties can also be evaluated 
exactly in the absence of free vortices [|l8| . For nonequi- 
librium situations where free vortices and antivortices are 
present, e.g. after a quench into the KT phase from a dis- 
ordered (high-temperature) initial condition, one can ar- 
gue as follows. In the late stages of coarsening, when the 
remaining free vortices and antivortices are widely sepa- 
rated, with a typical spacing S,{t), the spin wave theory 
can be used on scales much larger than the microscopic 
scale a but smaller than £,{t). For example, the calcula- 
tion of the dynamics of a single, widely-separated vortex- 
antivortex pair would proceed as at T = 0, but using the 
temperature-dependent function ps{T) and r(r) that in- 
corporate the effect of thermal fluctuations on smaller 
length scales. [Strictly one should use a scale-dependent 
spin wave stiffness and kinetic coefficient, evaluated at 
scale £,{t), but we are interested in the limit ^(i) -^ 00 
where these functions can be replaced by their infinite- 
scale limits, ps{T) and F(T)]. 

In this approach the only modification of the T = 
results would be the replacement of Ps and F, in the ex- 
pression for the annihilation time of a vortex-antivortex 
pair, by their T-dependent generalizations, as follows. 
Let ri and V2 be the positions of the vortex and antivor- 
tex, and r = r2 — ri be their relative separation. The 
equations of motion for ri and r2, at T = 0, read 

dvi 2ttPs^ 



, .dV2 



2'!rps 



(36) 



where "f(r) ~ (tt/F) ln(r/a) is the vortex mobility. Sub- 
tracting these gives the equation for the relative separa- 
tion r: 



^(r 



'di 



4:TTPs^ 

r. 



(37) 



Notice that the force —Airps/r corresponds to an effective 
potential Vcs = 2T^(r), where V{r) = 27rps ln(r/a) is the 
energy of a vortex-antivortex pair. 

How do we generalize this equation to T > 0? There 
are two distinct effects. Firstly, as we have discussed, 
thermally induced bound vortex-antivortex pairs renor- 
malize the large-scale spin-wave stiffness and kinetic co- 
efficient to temperature-dependent functions Ps{T) and 
F(T), with non-zero values Ps{Tkt) and T{Tkt) at the 
transition, i.e. the spin-wave theory describes the large- 
distance, large-time behavior at and below Tkt- Above 
Tkt, the spin wave theory breaks down due to the ther- 
mal nucleation of free vortices and antivortices. 

If the renormalization of ps and F were the only ef- 
fects of thermal fluctuations, we would conclude imme- 
diately that deterministic annihilation of a single vortex- 
antivortex pair would occur, as at T = 0, on a timescale 
t ~ Tq ln(ro/a), where tq is the initial separation. A sec- 
ond consequence of thermal noise, however, is diffusion 
of the vortex and antivortex. This means we have to add 
Langevin noise terms to (|3^). We anticipate that the 
strength of the noise will be proportional to 1/7 (r) on 
account of the fluctuation-dissipation theorem (Einstein 
relation) . Therefore we write a stochastic version of (|37| ) 
in the form 



rfr 



where 



l{r)- = -VU+^^)m, 



(6(i)e,(i')>-4i^'5y'5(i-0, 



(38) 



(39) 



and U{r) is a central potential which, it turns out, 
will differ by terms of order D from the effective po- 
tential Vcfi (r) — i-Kps In r suggested by (^ . We will 
determine U{r) from the condition that the stationary 
distribution Q{r) derived from ( p8[) satisfles Q{r) oc 
27rr exp[— 14ff(r)/2D]. The strength of the noise in ( |3q ) 
is 4D, rather than the usual 2D (where D — ksT) since 
the noise acting on r is the sum of independent noises of 
strength 2D acting on ri and r2. 

The noise term in (|38| ) is multiplicative noise in the 
Stratonovich sense. To determine the stationary distri- 
bution, and hence infer the form of U{r), we first write 
(Bq) in the canonical form 



|=A(r)+.9(r)|(t), 



(40) 



where 



1 

1, 



9 

7 



The corresponding Fokker-Planck equation is 
V ■[-AP + 2DgV{gP)]. 



dP 
'dt 



(41) 
(42) 

(43) 



Noting that A(r) is in the radial direction, i.e. A(r) = 
—g^{r) {drU)r, this equation can be reduced to a one- 
dimensional equation for the angle-averaged radial dis- 
tribution function, Q{r,t) — t"^^^ J dnP{r,i^,t), just as 



m section 



HA 



where 



dtQ = dr[g\drW)Q + 2Dgdr{gQ)], 



W{r) ^ U{r)~2D{d-l)lnr. 



(44) 



(45) 



Equation (W4j) can be recast as the one-dimensional 
Langevin equation 



dr 9, , dW , s ^, s. 

^,--9\r)^+g{r)m, 



(46) 



with {£,{t)£,{t')) = AD6{t - 1'), which in turn can be writ- 
ten as a Langevin equation with additive noise, via the 
change of variable dy — dr / g{r): 



dy 
dt 



dW 
dy 



m- 



(47) 



We can now determine the potential U{r) that 
corresponds to a 'physical' potential V{r). Since 
the effective potential for the dynamics of y is 
W{y), the stationary distribution for y is Py{y) oc 
exp[— W^(j/)/2Z?], implying that the radial distribu- 
tion function, Q{r) = Py{y){dy / dr), is given by 
Q{r) ex [g{r)Y^cyip[~W{r)/2D] which, using (H), im- 
plies Q{r) oc r''-iexp{-[C/(r) - 2i:>lng(r)]/2i:>}. The 
prefactor r''"^ is just the phase-space factor for the d- 
dimensional space: the required stationary distribution 
function for motion in a potential Vcsir) [—2V{r)\ is 
Q{r) ex r"^"^ exp[— V;3ff(r)/2D]. Equating these two re- 
sults for Q{r), and using 14ff(r) = 2V{r), where V{r) is 
the vortex- antivortex potential energy, determines U{r): 



U{r) = 2V{r) + 2D\ng{r). 
The final equation for y{t) becomes 



(48) 



dy ( dV 2Ddg 2{d-l)D\ dr ^, , ,^^, 
dt \ dr g dr r / dy 

The final step is to apply this (rather general) result 
to the 2D XY model. Inserting 7(r) = 70 ln(r/a) and 



m 



V{r) = aln(r/a), where 70 = tt/F and a = 27rp„, 
(p^), using dr/dy = g{r) — \/ \J^{r), and setting d — 2 
gives 



dt 



VW) 



D~a 



7S'^«"- 



(50) 



Consider now the size of the second term in the large 
bracket. Using ( [4l| ) gives 



1 dg 
9 dr 



1 din 7 

2 dr 



1 1 

2rln(r/a) r 



for r — > cx). This means that this term can be dropped 
in the calculation of the large-distance behavior of the 
system. Furthermore, integrating the relation dy/dr = 
\/j{r) gives, to leading order for large r, y = ry^j{r) — 
r-y/7^)hi(r/a), with corrections which are again of rela- 
tive order l/ln(r/a). To leading order, therefore, (|& 
reduces to (inserting a = 27: ps) 



dy 2D - A-Kps 



dt 



y 



-m, 



(51) 



with (C(i)C(t')> =4:Dd{t-t'). 

Apart from a doubling of the noise strength, D — > 2D, 
associated with the fact that we are dealing here with a 
two-body problem, the result (B^) is identical to equa- 
tion (g) derived in section II A for a particle moving in a 
logarithmic potential, except that the coordinate y corre- 
sponds to r-\/7oln(r/a) in the present context as a result 
of the scale-dependent mobility. Vortex-antivortex anni- 
hilation can be deemed to have occurred when r — a, 
the vortex core size, which corresponds to y = 0, so our 
previous results on persistence in a logarithmic potential 
can be applied directly to the persistence problem for 
a vortex-antivortex pair. In particular the parameter b, 
which controls the first-passage time distribution, Pi{t), 
and persistence probability, Po(i), according to dl4) and 
( llSl ) respectively, is obtained by setting A = -iirps , d = 2 
and Z? -> 2L> = 2kBT in (|) to give 



^_ 27Tp,iT) 

ksT 



(52) 



The ratio ps{T)/kBT is a decreasing function of T, 
and approaches the universal limit 2/tt for T — > Tkt- It 
follows that 6 is a decreasing function of T, diverging to 
infinity for T — > 0, and approaching the limiting value 
3 for T -^ Tkt- Recall that the persistence exponent is 
9 = (1 + 6)/2, so 9 varies continuously between 9{Q) = oo 
and 9 (Tkt) = 2, where the infinite limiting value at 
T = simply reflects the deterministic collapse in the 
absence of thermal noise. The persistence exponent for 
vortex-antivortex annihilation is therefore non-universal, 
depending continuously on T: 



keT ' 



(53) 



for T < Tkt- From (|4|) we see that the mean annihi- 
lation time is finite for all T < Tkt (since b > 3), and 
the characteristic annihilation time scales as t ~ i/p ^ 
rQln(ro/a) where vq is the initial separation. We con- 
clude that thermal fluctuations do not change the funda- 
mental relation between lengthscales and timescales de- 
duced earlier from zero-temperature considerations, and 
that the coarsening growth law, ^{t) ~ (t/lnt)^/^, holds 
for aU T < Tkt- 



D. Comparison with Simulation Data 

Recent simulation data |16|] support the conclusions of 
the above analysis. Since the method and results have 
already been presented in a short paper |lq| , we will 
just briefly recall the salient features here. The nearest- 



neighbor lattice Hamiltonian H 



E 



i.j) 



COS 



-9,) 



was simulated on square lattices of size L x L, for 12 < 
L < 48. Both uniform (all spins parallel) and disor- 
dered (randomly orientated spins) initial conditions were 
employed. In both cases, the 'time-dependent Binder cu- 
mulant' 1211,^1 



g{L,t) = 2- 






(54) 



where M{t) is the total magnetization at time t, was 
measured, and a finite-size scaling, g{L,t) = f[^{t)/L\, 
attempted, where ^(t) is a characteristic length scale at 
time t. 

In [|l6| , all simulations were at the KT temperature 
T = 0.90 Q. In earlier work, Luo et al. (24) have 
additionally studied some lower temperature within the 
KT phase. For uniform initial conditions, a good scal- 
ing collapse was obtained using a characteristic length 
scale S,{t) ~ t^/'^, with an exponent z that tends to 2 
(from slightly smaller values) for large L in accordance 
with the predictions of spin- wave theory. This is entirely 
reasonable since the uniform initial conditions contain 
no free vortices and subsequent thermal fluctuations can 
only produce bound vortex-antivortex pairs. 

For disordered initial conditions, a naive collapse us- 
ing ^(i) ~ t^'^- gives unreasonably large values of z, in 
the range 2.3 — 2.4, ||2j,|l^, for all temperatures stud- 
ied. Disordered initial conditions seed the system with 
free vortices and antivortices, so the analysis presented 
above suggests that a characteristic length scale ^(i) ~ 
(t/lni)^/^ is appropriate for this case. This form indeed 
yields an excellent data collapse. In practice, the form 
S,{t) = [t/\Ti(t/tQ)Y/'^ is used, with io a fitting parameter 
whose value is of order unity ||l^ . 

The data can, however, be presented in another way. 
The logic of finite-size dynamical scaling requires that the 
scaling variable can be written as t/r^L), where t(L) is 
the relaxation time of the finite-size system. For uniform 
initial conditions, the spin-wave theory gives r(i) ~ i^. 




which is equivalent to the previous scaling form, but for 
disordered initial conditions t(L) is given by the time for 
the last free vortex-antivortex pair, which are typically 
separated by a distance of order L, to annihilate, i.e. 
t{L) ~ L^ InL. A scaling collapse using i/[L^ ln(L/Lo)] 
as scaling variable is presented in Figure 1, with Lq = lA 
where Lq is a short distance cut-off expected (and found) 
to be of order unity. The scaling collapse is excellent. For 
large enough L and t, this method of plotting the data 
and the previous method, using t/[L^ln(i/io)] as scaling 
variable, will be indistinguishable. For finite times and 
sizes, however, they differ slightly, and the new method 
gives an improved collapse. 
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0.6 A 
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FIG. 1. Scaling plot for the time-dependent Binder param- 
eter, starting from a disordered initial condition, for system 
sizes L = 12, 16, 24, 32, 48 and Lo = lA 



E. Discussion 

The fact that £,{t) depends on the initial conditions is 
surprising from the viewpoint of conventional nonequi- 
librium critical dynamics [Q, according to which ^(i) ~ 
t^'^, where z is the dynamical exponent for equilibrium 
critical dynamics. An equivalent statement in the con- 
text of finite-size dynamical scaling is that the relaxation 
time grows with L as t{L) ^ L^, independent of ini- 
tial conditions. While this seems to be true for most 
phase transitions (e.g. the 2D Ising model [g5|) we hope 
we have convinced the reader that it docs not hold for 
the 2D XY model, at and below Tkt, for disordered ini- 
tial conditions. We conjecture that the breakdown of the 
standard field-theoretic methods ||ll[] in this case is due 
to the key role played by vortex configurations which, 
due to their non-perturbative character, are not accessi- 
ble to perturbative methods based on a 4 — e expansion. 
The latter involves perturbing around a Gaussian theory, 
which cannot support topological defects. It would be in- 



teresting to investigate whether similar resuhs obtain in 
other defect-driven phase transitions. 
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